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Abstract 


This  work  demonstrates  the  first  phases  of  a  newly  proposed  gun  propellant  formulation 
process  that  will  minimize  life-cycle  costs  through  science-based  design.  This  new  approach 
proposes  maximal  use  of  modeling  and  simulation  in  the  earliest  phases  of  the  developmental 
cycle  to  screen  candidate  formulations,  resulting  in  elimination  of  probable  poor  performers  and 
identification  of  the  most  promising  test  candidates.  The  screening  and  identification  of 
propellant  formulations  are  demonstrated  under  the  assumption  of  a  specific  weapon  platform 
and  user  requirements.  The  process  of  selecting  a  propellant  for  the  assumed  gun  system 
application  has  been  distilled  into  measurable  steps,  leading  from  a  set  of  candidate  materials, 
through  logical  and  numerical  filters,  to  a  shorter  list  of  energetic  materials  demonstrated  as 
viable  weapon  platform  choices.  Environmental  filtering  and  performance  modeling  are  used 
to  screen  propellants  through  a  well-defined  sequence  of  tests  designed  to  weed  out  materials  not 
meeting  safety,  energy,  or  manufacturability  standards.  Because  much  of  the  testing  is  performed 
by  computer  modeling,  the  gun  systems  and  energetic  materials  need  not  be  present  (or  even 
existent)  in  order  to  be  described  and  matched  against  performance  requirements  for  future 
applications.  The  calculations  demonstrate  that  utilizing  computer  models  rather  than  physical 
testing  in  the  early  developmental  stages  of  the  formulation  process  can  produce  enormous 
savings  in  labor,  material,  and  environmental  costs,  along  with  a  tremendous  reduction  in  the 
time  required  to  select  a  “best  candidate”  propellant. 
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1.  Introduction 


The  attainment  of  proposed  future  war-fighting  capabilities  [1]  wiU  depend  on  the 
performance  of  new  weapons  platforms  and  the  energetic  materials  that  are  used  by  them.  In 
order  to  efficiently  and  effectively  meet  the  requirements  of  the  future  weapons  systems,  a 
mechanism  must  be  in  place  for  the  rapid  design  and  manufacture  of  advanced  chemical  fuels. 
Presently,  there  is  no  single  formal  procedure  in  place  for  the  design  and  development  of  a  new 
propellant;  rather,  the  overall  procedure  includes  any  number  of  steps  leading  to  deployment  of  a 
new  energetic  material.  The  inclusion  of  steps  varies  among  individuals  according  to 
experience,  philosophy,  and  training.  Also,  the  overall  procedure  relies  heavily  on 
experimentation  and  measurement  and  often  proceeds  through  lengthy  cut-and-try  processes. 
Significant  waste  streams  in  terms  of  time,  manpower,  and  materials  are  generated  from 
implementation  of  such  unstructured  developmental  procedures,  and  they  often  exact 
unacceptably  large  time  and  environmental  penalties.  Current  fiscal  realities  and  increased 
environmental  constraints  necessitate  a  revision  in  existing  strategies  for  energetic  material 
design,  development,  and  manufacture.  The  most  obvious  revision  that  must  be  made  is  to 
minimize  the  waste  created  during  the  process.  In  order  to  do  this,  a  formal  mechanism  must  be 
established  to  estimate  the  waste  associated  with  each  part  of  the  procedure  and,  thus,  assign 
costs  to  each  step  in  the  mechanism.  The  assignment  of  costs  will  depend  upon  the  priorities  that 
drive  the  development  of  the  material  and  could  include  pecuniary,  environmental,  or 
performance  considerations.  Then  refinements  to  the  steps  and  overall  mechanism  can  be  made 
to  minimize  the  cost  and  maximize  the  efficiency  of  the  process. 

The  Strategic  Environmental  Research  and  Development  Program  (SERDP),  in  recognizing 
the  need  to  counter  extreme  waste  generation  in  energetic  materials  development,  is  supporting  a 
project  whose  goal  is  to  establish  a  new  developmental  procedure  that  will  minimize  waste 
generation  throughout  the  entire  process.  The  project  is  currently  limited  to  gun  propellant 
formulations;  however,  the  same  developmental  procedures  could  be  applied  to  other  classes  of 
energetic  materials,  including  rocket  propellants,  explosives,  or  pyrotechnics.  The 
developmental  procedure  that  is  being  advanced  under  this  project  is  model-based.  In  other 
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words,  the  goals  of  the  project  can  ideally  be  achieved  through  replacement  (where  possible)  of 
measurement  and  testing  with  computer  modeling  or  through  optimal  reduction  of 
experimentation  by  the  exercise  of  predictive  methodologies.  This  new  approach  to  gun 
propellant  formulation  is  outlined  in  a  parallel  report  [2]  that  describes  a  science-based  activity 
where  the  earliest  possible  steps  utilize  modeling  and  simulation  of  all  stages  that  might  occur 
during  development  of  a  gun  propellant.  Although  modeling  of  the  complete  process  cannot  be 
implemented  at  this  time  since  there  are  several  areas  in  which  models  do  not  currently  exist,  the 
work  described  here  is  an  attempt  to  exercise  the  developmental  procedure  using  available 
models. 

Part  of  the  propellant  assessment  involves  determination  of  whether  its  performance  in  a 
cannon  will  meet  the  military  requirements  of  projectile  velocity  (which  governs  effective  range 
of  the  weapon)  within  the  limitations  of  pressure  and  acceleration  peculiar  to  the  gun  system  and 
projectile.  Predictive  performance  models  [3—5]  have  been  available  in  the  military  for  many 
years  and  are  well-suited  for  computing  expected  solutions  of  the  interior  ballistic  cycle  for 
known  (measurable)  guns  and  for  gun  systems  as  yet  only  imagined.  These  computer  models 
require  a  priori  information  describing  the  physical  parameters  of  the  gun  chamber  and  tube,  of 
the  projectile,  of  the  propellant  charge,  and  of  critical  interactions  during  the  gun  firing.  This 
report  explains  how  to  use  available  computer  programs  to  model  the  performance  testing  of  a 
proposed  set  of  gun  propellants  in  a  large-caliber  artillery  weapon.  Additionally,  since  the 
proposed  developmental  approach  incorporates  environmental  screening  procedures  at  each 
stage  of  the  development  process,  we  have  attempted  to  include  such  in  our  demonstration.  A 
prototypical  simulation  was  performed  in  which  minimal  user  requirements  for  maximum 
pressure  and  minimum  muzzle  energy  were  assigned  for  a  specific  weapon  platform  using  a  set 
of  10  existing  propellant  formulations.  Using  performance  modeling  and  database  information 
on  environmental  hazards,  6  of  the  10  formulations  were  eliminated  as  candidates  due  to 
environmental  constraints  and  failure  to  meet  energy  requirements.  This  demonstration  shows 
that  application  of  computer  models  based  on  scientific  principles  during  gun  propellant  design 
and  formulation  will  result  in  cost  savings  through  decreases  in  labor,  equipment,  and 
environmental  issues. 
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2.  Performance  Testing  Description 


2.1  Propellant  Database.  A  generalized  flow  chart  of  the  necessary  steps  in  testing  the 
performance  of  a  new  gun  propellant  is  shown  in  Figure  1,  representing  the  early  screening  steps 
of  the  science-based  design  process  outlined  in  the  study  of  Miller,  Rice,  and  Cramer  [2].  The 
Master  Optimization  Algorithm  (MOPTA)  starts  with  a  database  of  new  and/or  known 
propellants  to  be  tested  and  screened  for  suitability  to  a  “virtual  gun  system.”  This  gun  system 
may  be  a  computer  model  of  a  known  weapon  or  a  model  of  a  completely  theoretical  gun  system 
to  be  analyzed  for  possible  manufacture  and  deployment.  From  its  chemical  and  thermodynamic 
description,  each  propellant  is  to  be  characterized  in  the  weapon  without  ever  having  physically 
been  placed  in  the  chamber,  or  perhaps  without  ever  having  been  actually  created  in  a  laboratory. 
The  portion  of  the  algorithm  to  be  described  in  this  report  is  highlighted  with  a  gray  background 
in  Figure  1 — ^the  other  blocks  contain  steps  where  either  models  are  not  yet  available  or  else  the 
authors  do  not  have  the  expertise  necessary  to  use  current  programs  or  databases. 

2.2  Ingredient  Environmental  Screening.  Those  propellants  containing  known 
environmental  hazards  may  be  screened  out  of  the  performance  testing  process  here.  Legacy 
propellants  may  include  ingredients  that  have  been  identified,  in  the  years  since  original 
production,  as  carcinogens  or  as  substances  toxic  to  one  or  more  environmental  systems.  This 
could  include  not  only  the  ingredients  themselves  but  also  the  manufacturing  processes  involved 
in  creating  the  material. 

23  Thermodynamic  Properties.  Each  propellant  must  be  described  by  a  suite  of 
measurable  physical  parameters  in  order  to  be  acceptable  to  the  interior  ballistic  (EB)  computer 
programs.  Required  values  include  isochoric  flame  temperature,  density,  internal  chemical 
energy,  and  others  depending  on  which  computer  code  will  be  used  for  performance  modeling. 
Several  programs  are  available  to  estimate  these  parameters  (except  for  density  and  burning  rate) 
from  the  initial  chemical  description.  These  programs  include  CHEETAH  [6],  BLAKE  [7],  and 
TRAN72  [8]  and  CEA  [9].  Part  of  the  output  generated  by  these  codes  can  include  a  list  of 
combustion  products  (both  gaseous  and  condensed),  which  can  be  scanned  to  determine 
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DVANCED  EVALUATION  MODELS 


Figure  1.  Master  Optimization  Algorithm. 
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if  toxic  or  corrosive  materials  are  predicted.  Propellant  density  and  burning  rates  must  be 
physically  measured  or  estimated  from  other  sources. 

2.4  Estimate  of  Available  Energy.  Given  a  gun  system  model  and  the  thermodynamic 
properties  of  each  propellant,  an  estimate  can  be  made  of  how  much  energetic  material  can  be 
placed  in  the  gun  chamber.  A  calculation  can  then  be  made  to  determine  if  enough  kinetic 
energy  can  be  transferred  to  the  projectile  in  order  to  meet  the  exit  velocity  requirements.  By 
assuming  that  the  highest  allowable  pressure  is  continuously  maintamed  by  propellant 
combustion  (until  the  charge  is  exhausted),  an  estimate  of  maximum  achievable  projectile  exit 
velocity  can  be  calculated.  If  this  velocity  does  not  meet  or  exceed  user  requirements,  then  the 
propellant  will  not  be  a  candidate  for  further  testing.  Computer  codes  that  can  estimate  this 
maximum  are  IBHVG2  [10]  and  CONPRESS  [1 1]. 

2.5  Propellant  Grain  Geometry.  Actual  energy  transferred  to  a  projectile  will  be 
somewhat  less  than  that  calculated  in  the  previous  step  because  propellant  grain  design  caimot 
provide  a  continuous,  unwavering  pressure  level  for  the  life  of  the  burning  charge.  Standard 
grain  geometries  (19-perf,  7-perf,  and  single-perf,  shown  in  Figure  2)  can  be  parametrically 
computed  and  tested  within  the  gun  system  model.  In  general,  a  multiperforated  grain  design  is 
more  efficient  in  transferring  energy  to  a  projectile;  this  rule  of  thumb  may  be  untrae  for  certain 
“virtual”  gun  systems  where  databases  may  be  varied  to  experiment  with  nontraditional 
dynamics  in  chamber  sizes,  projectile/bore  interactions,  and  other  physical  parameters.  Interior 
ballistic  calculations  can  determine  if  one  or  more  propellant  grain  designs  used  as  charges  in  the 
real  (or  future)  gun  system  can  meet  requirements  for  muzzle  energy  and  complete  charge 
burnout  within  maximum  pressure  limitations.  A  first-pass  IB  program  is  IBHVG2;  if 
predictions  of  pressure  waves  within  the  gun  chamber  are  required,  then  XNOVAKTC  [12]  or 
NGEN  [13]  may  be  used. 

2.6  Low-Zone  Charge  Applications.  Artillery  guns  are  routinely  loaded  with  a 
less-than-maximum  charge  when  targets  are  located  at  relatively  short  ranges.  These  “zoned” 
applications  allow  optimum  trajectories  for  higher  kill  probabilities.  Lowest-zoned  charges  are 
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19Perf. 


7Perf. 


1  Perf. 


Figure  2.  Propellant  Grain  Types. 

still  required  to  bum  completely  before  the  projectile  exits  the  gun  muzzle;  unbumed  and 
partially  burned  propellant  particles  and  gases  can  be  major  causes  of  secondary  muzzle  blast 
and  flash.  If  low-zone  applications  are  defined,  then  the  grain  geometry  meeting  minimum 
energy  requirements  at  maximum  charge  weights  must  still  show  complete  burnout  when  used  at 
lower  loading  densities.  "Die  same  IB  programs  that  previously  calculated  acceptable  propellant 
grain  geometries  may  also  be  used  here. 

2.7  Propellant  Suitability.  Other  conditions  may  screen  out  propellants  that  meet  muzzle 
energy  requirements.  Cost,  manufacturing  limitations,  burning  temperature  (a  major  factor  in 
gun  tube  life  expectancy),  grain  size,  or  any  number  of  user-defined  stipulations  may  preclude 
propellants  from  final  consideration.  More  likely,  the  list  of  remaining  propellants  will  be 
ranked  according  to  importance  of  these  additional  factors. 

3.  Example  Performance  Screening 

3.1  User  Requirements.  As  a  test  system,  the  Navy  5-inch  54-caliber  gun  was  chosen  to  be 
the  vehicle  for  exercising  the  propellant  performance  analysis.  User-supplied  guidelines 
included  a  maximum  pressure  limitation  of  65  ksi  (448  MPa)  and  a  minimum  muzzle  energy  of 
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18  MJ  with  a  1 10-lb  (49.9  kg)  projectile.  From  the  kinetic  energy  formula  (KE  =  V2  MV^),  it  can 
be  determined  that  projectile  velocity  needs  to  be  at  least  848  m/s  to  meet  the  requirement.  A 
general  description  of  the  weapon  specifies  a  chamber  volume  of  approximately  0.0151  m 
(919  in^)  and  a  projectile  travel  of  6.84  m  (269  inches). 

3.2  Propellant  Database.  A  mixture  of  old  and  new  propellant  formulations,  from  both 
Army  and  Navy  sources,  provided  a  well-rounded  sample  for  the  analysis.  All  propellants  are 
real — ^no  theoretical  materials  have  been  added  to  the  study,  although  it  certainly  could  be  done 
in  further  tests.  The  standard  Army  propellants  include  Ml,  MIO,  M26E1,  M30A1,  M31E1,  and 
JA2;  Navy  propellants  are  NACO,  EX99  (called  LOVA  in  the  reference  noted  in  Table  1), 
HELP42,  and  BAMO/AMMO  (abbreviated  as  BAMO  in  the  rest  of  diis  report).  No  particular 
method  was  used  for  propellant  selection  other  than  an  attempt  to  include  the  several  classes 
(single-base,  double-base,  and  triple-base  propellants),  along  with  some  newer  formulations  not 
fitting  neatly  into  such  descriptions. 

3.3  Ingredient  Environmental  Screening.  Table  1  shows  the  proposed  propellants  and 
fists  their  ingredients.  A  quick  scan  reveals  that  NACO  contains  lead  carbonate  iPb2C04) — a 
known  toxic  substance  and  suspected  carcinogen  [14].  For  this  reason,  NACO  was  dropped 
fi-om  further  analysis.  (If  a  cost  of  manufacture  with  suitable  biohazard  controls  and  increased 
demilitarization  expenses  would  be  computed,  then  such  a  material  could  be  included  in  the 
complete  performance  appraisal;  the  extra  factors  would  be  added  to  a  final  figure  for  overall  fife 
cycle  costs.) 

Dinitrotoluene  (DMT)  is  an  ingredient  of  the  Army  Ml  propellant.  The  fiitemational 
Chemical  Safety  Card  (ICSC  0727)  [15]  describes  this  substance  as  able  to  be  absorbed  through 
human  skin  and  ingested  through  inhalation  of  fumes  when  it  is  heated.  For  humans,  it  is  an 
irritant  to  the  eyes  and  skin  and  may  cause  effects  on  the  central  nervous  system,  cardiovascular 
system,  and  blood.  It  is  also  extremely  toxic  to  aquatic  organisms.  The  environmental  costs 
from  DNT  can  be  very  high;  as  an  example  of  the  type  of  chemical  ingredient  to  be  avoided, 
mixtures  with  DNT  (specifically  Ml)  will  be  removed  from  the  fist  of  propellants  under 
consideration. 
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Table  1.  Ingredients  of  Propellants 


Ingredients 

Reference 

Ml 

NC1315  DMT  DBP  DPA  H20  ALC 

Freedman  [16],  p.  121 

MIO 

NC1315  DPA  KS  H20  ALC 

Freedman  [16],  p.  121 

M26E1 

NC1315  NG  EC  H26  ALC  C 

Freedman  [16],  p.  122 

M30A1 

NC1260  NG  NQ  EC  KS  ALC  C 

M31E1 

NC1260  NG  NQ  DBP  NDPA  KS  ALC  C 

Freedman  [16],  p.  122 

NACO 

NC1200  BS  EC  KS  PB2C04  H20  ALC 

Freedman  [16],  p.  122 

JA2 

NC1304  NG  DEGDN  AKAR2  MGO  H20  C 

Miller  [17] 

BAMO 

RDX  BAMO  AMMO 

Almeyda  [18] 

HELP-42 

HZTZ  RDX  NCI 260  BDNPF  BDNPA  EC 

Cramer  [19] 

EX-99 

RDX  NCI 260  CAB  BDNPF  BDNPA  EC 

Cramer  [19] 

The  ingredient  definitions  in  Table  1  and  in  text  are  as  follows: 


AKAR2 

Akardite  E 

ALC 

Ethanol  (ethyl  alcohol) 

AMMO 

Azidomethylenemethyl  oxetane 

BAMO 

Azidomethyl  oxetane 

BDNPA 

Bisdinitropropyl  acetal 

BDNPF 

Bisdinitropropyl  formal 

BS 

Barium  sulfate 

C 

Graphite 

CAB 

Cellulose  acetate  butyrate 

DBP 

Butyl  phthalate 

DEGDN 

Diethyleneglycol  dinitrate 

DNT 

Dinotrotoluene 

DPA 

Diphenylamine 

EC 

Ethyl  centralite 

H20 

Water 

HZTZ 

Bisdihydrazinotetrazine 

KS 

Potassium  sulfate 

MGO 

Magnesium  oxide 

NCxx.yy 

Nitrocellulose  with  a  nitration  percentage  of  xx.yy% 

NDPA 

2-nitrodiphenylamine 

NG 

Glyceryl  trinitrate 

NQ 

Nitroguanidine 

PB2C04 

Basic  lead  oxide 

RDX 

Cyclo- 1,3,5  -trimethylene-2, 4,6-trinitramine 
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3.4  Thermodynamic  Properties.  The  U.S.  Army  Research,  Development,  and  Engineering 
Center  (ARDEC)  Report  “Interior  Ballistics  Firing  Data  Library  of  Closed  Breech,  Single 
Combustion  Chamber  Tank  Guns,  Artillery  Guns  and  Howitzers  Cannon,”  by  author  Frank  J. 
Virginia,  Jr.  [20],  is  a  valuable  source  of  existing  military  propellant  data  including  many 
measured  burning  rates.  All  data  are  taken  from  testing  of  individual  propellant  lots.  A  second 
reference  is  the  chapter  “Thermodynamic  Properties  of  Military  Gun  Propellants”  written  by  E. 
Freedman  in  the  book  Gun  Propulsion  Technology  [16],  from  which  many  of  the  ingredient 
listings  in  Table  1  have  been  taken.  The  data  in  this  reference  pertain  to  generic  compositions 
rather  than  the  variations  found  in  individual  lots  of  propellant. 

Table  2  contains  thermodynamic  values  required  to  model  the  remaining  propellants.  As 
shown,  there  is  a  wide  range  of  isochoric  flame  temperatures  (from  lowest  at  2,555  K  to  highest 
at  3,410  K).  Estimated  temperature  within  the  gun  chamber  will  affect  modeling  of  product  gas 
expansion  due  to  the  (simplified)  nonideal  equation  of  state  p(V-b)=RT  [21],  where  T  is  a  direct 
factor  in  the  pressure  calculation;  in  turn,  the  force  of  pressime  on  the  base  of  the  projectile  will 
govern  acceleration  through  F=ma.  Propellant  burning  rates  will  also  be  a  major  factor;  rates  in 
this  study  were  gained  either  from  the  aforementioned  Virginia  reference  or  from  other  sources 
as  noted  in  Almeyda  [18]  and  Cramer  [19].  The  propellant  property  labeled  as  Force  in  Table  2 
is  a  measure  of  the  energy  contained  in  a  unit  weight  of  material. 


Table  2.  Thermodynamic  Properties  of  Propellants 


Propellant 

Density 

(g/cm^) 

Gamma 

Temperature 

(K) 

Covolume 

(cm^/g) 

Force 

(J/g) 

MIO 

1.669 

1.2342 

3,000 

1.0025 

1013.3 

M26E1 

1.611 

1.2384 

3,132 

1.0383 

1085.0 

M30A1 

1.683 

1.2375 

3,036 

1.0524 

1073.4 

M31E1 

1.642 

1.2580 

2,574 

1.1048 

973.8 

JA2 

1.661 

1.2255 

3,410 

0.9978 

1130.2 

BAMO 

1.639 

1.2738 

2,863 

1.2218 

1191.3 

HELP-42 

1.600 

1.2780 

2,555 

1.22 

1063.0 

EX-99 

1.660 

1.2900 

3,019 

1.122 

1159.0 
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If  there  were  no  published  or  experimental  thermodynamic  values,  computer  programs  are 
available  to  calculate  estimates.  CHEETAH,  BLAKE,  and  TRAN72  have  all  been  used  for  this 
purpose.  As  an  example,  the  formulation  for  Ml  (extracted  from  the  Freedman  chapter,  and 
shown  in  Table  3)  has  been  used  for  an  input  to  the  BLAKE  program.  The  output  is  reproduced 
in  Appendix  A;  the  thermodynamic  values  are  listed  in  Table  4,  along  with  comparable  data  from 
the  Virginia  report.  The  fluctuation  of  values  is  not  unexpected  because  any  particular  lot  of  Ml 
may  vary  slightly  from  the  composition  of  the  generic  propellant — ^manufacturing  methods  are 
not  so  precise  as  to  guarantee  final  composition  in  a  production  lot  of  propellant  to  within  a 
fraction  of  a  percent. 


Table  3.  Ingredients  and  Their  Relative  Percentages  for  Ml  Propellant 


Ingredient 

BLAKE  Abbreviation 

Nitrocellulose 
(13.15%  Nitration) 

NC1315 

83.11 

Dinitrotoluene 

DNT 

9.77 

Butyl  Phthalate 

DBP 

4.89 

Diphenylamine 

DPA 

0.98 

Water 

H20 

0.50 

Ethanol 

ALC 

0.75 

Table  4.  Ml  Thermodynamic  Data  Calculated  by  BLAKE  and  From  Virginia  Reference 


Source 

Lot  Number 

Gamma 

Temperature 

(K) 

Covolumc 

(cm^/g) 

Blake 

(Generic) 

1.2675 

2,447 

1.105 

919.2 

Virginia 

JAN-P-309 

1.2593 

2,417 

1.1044 

911.7 

Part  of  the  calculation  of  thermodynamic  parameters  includes  estimation  of  final  combustion 
gases  and  their  relative  amounts.  A  section  of  Appendix  A  titled  “Constituent  Concentrations” 
shows  gases  (and  possible  solid  or  liquid  products)  expected  to  be  created  by  each  mole  of 
burned  energetic  material.  If  relatively  large  amounts  of  toxic  or  corrosive  products  are  noted, 
the  propellant  could  be  dropped  from  further  analysis  or  else  flagged  for  special  consideration 
later  in  the  study. 
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The  computation  of  burning  rate  parameters  is  not  yet  so  straightforward  as  is  the  case  with 
thermodynamic  values,  since  there  is  no  computer  program  to  calculate  rates  from  chemical 
composition.  Obviously  the  best  source  is  through  experimental  data  from  closed-bomb  analysis 
of  the  actual  propellant  [22-24].  Estimates  of  burning  rates  and  propellant  densities  could  be 
taken  from  a  member  of  the  “family”  of  propellants  if  the  formulation  of  the  unknown  material  is 
similar  (within  a  percentage  point  or  two)  to  the  major  ingredients  of  a  known  propellant,  but  the 
user  should  be  aware  that  it  is  only  an  estimate  and  needs  to  be  proven  through  experimental 
means. 


3.5  Estimate  of  Available  Energy.  IBHVG2  computations  were  used  to  evaluate  the 
potential  energy  of  each  propellant  in  the  gun  system  model.  An  assumption  made  for  this  step 
was  that  there  would  be  a  constant  charge  weight  for  the  various  energetic  materials.  Since  most 
fielded  granular  charges  have  maximum  loading  densities  (amount  of  propellant  mass  per  unit 
chamber  volume)  in  the  range  of  0.80  g/cm^  to  0.90  g/cm^,  a  suitable  value  for  this  study  was 
chosen  at  0.85  g/cm^  so  that  the  interior  ballistics  program  now  has  aU  the  necessary  parameters 
to  proceed  with  a  constant-pressure  calculation.  The  purpose  of  this  step  is  to  quickly  determine 
whether  the  charge  could  meet  the  minimum  requirement  of  projectile  exit  velocity  while 
keeping  chamber  pressure  no  higher  than  the  user-defined  maximum.  The  program  artificially 
ties  breech  pressure  to  a  given  value  (by  converting  the  necessary  amount  of  propellant  to  gas  at 
each  time  step)  until  the  charge  is  con^letely  burned;  then  the  gases  are  allowed  to  continue 
expanding  (and  accelerating  the  projectile)  until  maximum  travel  is  accomplished.  This  process 
approximates  what  would  be  the  perfect  combination  of  propellant  surface  area  and  burning  rate 
in  order  to  transfer  maximum  energy  to  the  projectile.  CONPRESS  does  not  require  estimates 
for  either  grain  geometry  or  burning  rate;  IBHVG2  requires  either  burning  rate  or  grain 
geometry,  although  this  is  just  a  formality — ^both  can  be  made  to  vary  in  order  to  complete  the 
calculation. 

Computed  projectile  exit  velocities  using  IBHVG2  are  listed  in  Table  5  for  each  of  the 
remaining  candidate  propellants.  The  user-required  velocity  is  848  m/s;  each  velocity  in  the 
table  is  compared  to  that  value  via  a  percentage  of  minimum  requirement.  All  considered 
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Table  5.  Constant-Pressure  Calculations  of  Exit  Velocity 


Propellant 

Exit  Velocity 
(m/s) 

Comparison 

(%) 

MIO 

867 

102.2 

M26E1 

858 

101.2 

M30A1 

886 

104.5 

M31E1 

825 

97.3 

JA2 

902 

106.4 

BAMO 

914 

107.8 

HELP-42 

879 

103.7 

EX-99 

894 

105.4 

propellants  except  M31E1  attained  over  100%  of  minimum  and  will  be  passed  on  to  the  next  step 
in  the  study. 

A  sample  output  from  the  constant-pressure  computation  for  M26E1  is  included  as 
Appendix  B  of  this  report. 

3.6  Propellant  Grain  Geometry.  A  parametric  feature  allows  IBHVG2  to  vary  propellant 
grain  dimensions  in  order  to  calculate  how  each  geometry  will  perform  in  the  modeled  gun 
system.  By  adjusting  one  grain  measurement  (in  this  case  the  web,  or  bum-through  distance 
between  grain  perforations  and  between  the  outer  grain  surface  and  closest  perforations)  the 
program  can  compute  entire  balhstic  cycles  for  a  series  of  grain  sizes  and  can  find  the  maximum 
pressure  and  the  expected  projectile  exit  velocity  for  each  situation.  Appendix  C  is  the  printout 
from  the  performance  calculations  for  M30A1;  it  includes  a  table  of  input  and  output  values  for 
each  set  of  dimensions.  A  short  summary  of  that  information  is  shown  in  Table  6,  where  the 
computed  maximum  pressure  is  closest  to  448  MPa  at  a  web  size  of  3.06  mm  for  the  grain  with 
19  perforations;  for  the  7-perf  grain  the  targeted  pressure  is  reached  at  a  web  size  of  3.38  mm; 
the  single-perf  simulations  show  maximum  pressure  of  448  MPa  at  the  web  thickness  of  4.72 
mm.  Projectile  exit  velocities  are  861  m/s,  856  m/s,  and  812  m/s,  respectively,  showing  that  the 
multiperforated  grains  transferred  more  energy  to  the  projectile  than  did  the  mono-perforated 
grain.  Both  the  19-perf  and  7-perf  solutions  exceeded  minimum  required  projectile  velocity. 
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while  the  single  did  not.  Other  grain  dimension  sizes  (length  18.3  mm  and  perforation  diameter 
0.457  mm)  were  held  constant  during  the  parametric  variations  of  web — ^this  means  the 
computed  grain  diameter  of  the  19-perf  solution  is 

3.06  *  6  +  0.457  *  5  =  20.6  mm, 

which  would  make  the  grain  length  slightly  shorter  than  its  width.  By  increasing  the  perforation 
diameter,  the  grain  could  be  longer  (perforation  length  is  restricted  by  the  ability  of  combustion 
gases  to  escape  confinement,  dependent  on  diameter  and  expected  burning  rates  of  the 
propellant) — ^but  this  is  an  exercise  best  left  to  expert  ballisticians.  For  current  purposes,  it  is 
enough  to  realize  that  the  multiperforated  M30A1  grains  can  satisfy  system  requirements. 

Table  6  also  contains  information  about  pressure  felt  on  the  projectile  base  at  the  time  of 
muzzle  exit.  The  “Z@Exit”  column  is  an  estimate  of  the  mass  fraction  burned  at  projectile  exit. 
The  last  column  is  an  estimated  length  of  projectile  travel  when  computed  charge  burnout 
occurs,  hi  the  bottom  two  rows  of  single-perf  data,  the  charge  is  not  expected  to  bum  completely 
by  the  time  the  projectile  exits  the  gun  muzzle;  this  could  be  a  clue  to  the  possibility  of 
unacceptable  secondary  muzzle  blast  and  flash  should  the  charge  consist  of  this  grain  geometry. 

For  the  19-perf  granular  solution  in  Table  6,  projectile  travel  is  only  slightly  more  than  2/3 
maximum  allowable — suggesting  that  there  is  yet  more  flexibility  in  charge  manipulation  in 
order  to  get  even  higher  exit  velocities  within  user  limitations.  For  example,  with  a  slightly 
larger  loading  density  (assuming  more  propellant  can  be  loaded  into  the  chamber)  and  a  larger 
grain  web  in  order  to  keep  down  the  maximum  pressure  for  the  increased  charge,  the  results  are 
as  in  Table  7.  As  charge  mass  increases  from  12.50  kg  to  13.50  kg,  the  web  size  is  increased 
from  3.0168  mm  to  3.3527  mm  in  order  to  keep  maximum  pressure  at  the  448-MPa  limitation, 
(hicreased  grain  web  size  results  in  a  lower  surface-to-mass  ratio;  thus  the  pressure  peak  can  be 
maintained  even  though  the  greater  charge  mass  transfers  more  energy  to  the 
projectile — assuming  all  the  propellant  still  bums  out  within  the  gun  tube.)  Loading  densities  for 
the  different  charge  weights  are  printed  in  the  second  column.  Values  in  the  right-most  column. 
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Table  6.  Computed  M30A1  Simulations  With  Different  Grain  Geometries  and  Web  Sizes 


3.00 

470.129 

868.56 

84.350 

1.000 

3.940 

3.02 

463.022 

866.16 

84.667 

1.000 

4.274 

3.04 

456.125 

863.60 

84.946 

1.000 

4.218 

3.06 

449.431 

860.93 

85.206 

1.000 

4.553 

mm 

442.933 

858.40 

85.503 

1.000 

4.495 

436.620 

855.71 

85.764 

1.000 

4.832 

3.12 

430.488 

853.11 

86.053 

1.000 

4.773 

3.14 

424.527 

850.57 

86.356 

1.000 

5.111 

3.16 

418.731 

847.84 

86.613 

1.000 

5.051 

3.18 

86.902 

1.000 

5.390 

3.20 

1.000 

5.328 

3.30 

458.389 

860.98 

85.192 

1.000 

4.700 

3.32 

453.193 

858.73 

85.494 

1.000 

5.050 

3.34 

448.116 

856.24 

85.729 

1.000 

4.999 

3.36 

443.157 

853.78 

85.979 

1.000 

5.351 

3.38 

438.312 

851.36 

86.244 

1.000 

5.298 

3.40 

433.576 

848.90 

86.498 

1.000 

5.651 

3.42 

428.946 

846.55 

86.795 

1.000 

5.598 

3.44 

424.419 

87.033 

1.000 

5.952 

3.46 

419.990 

841.55 

87.301 

1.000 

5.897 

3.48 

415.657 

839.03 

87.555 

1.000 

6.252 

836.59 

87.840 

1.000 

6.609 

4.60 

464.580 

r824.68 

89.090 

1.000 

5.530 

4.62 

461.798 

822.66 

89.303 

1.000 

5.498 

4.64 

459.051 

820.62 

89.473 

1.000 

5.860 

4.66 

456.340 

818.59 

89.694 

1.000 

5.828 

4.68 

453.663 

816.55 

1.000 

6.192 

4.70 

451.018 

814.50 

90.133 

1.000 

6.158 

4.72 

448.406 

812.44 

90.351 

1.000 

6.125 

4.74 

445.826 

810.37 

90.607 

1.000 

6.490 

4.76 

443.278 

808.29 

90.788 

1.000 

6.456 

4.78 

440.759 

806.21 

90.854 

0.999 

6.840 

4.80 

438.270 

804.13 

90.480 

0.995 

6.840 

Table  7.  M30A1  Charge  Weight  Variations  for  19-Perf  Granulation 


Charge 

Weight 

(kg) 

Loading 

Density 

(g/cm^) 

Web 

(mm) 

Maximum 

Pressure 

(MPa) 

Z 

Trav  @  B.O. 
(m) 

12.50 

0.838 

3.017 

448 

857.1 

1.0 

4.410 

12.75 

0.855 

3.097 

448 

862.4 

1.0 

4.473 

13.00 

0.872 

3.179 

448 

867.3 

1.0 

5.144 

13.25 

0.888 

3.264 

448 

871.6 

1.0 

5.561 

13.50 

0.905 

3.353 

448 

875.2 

1.0 

6.070 

projectile  travel  at  the  point  of  charge  burnout,  grow  larger  as  the  charge  weight  and  web  size 
increase — ^since  there  is  no  change  in  maximum  pressure  it  requires  a  longer  burning  cycle  to 
completely  consume  the  charge.  The  IBHVG2  computation  for  this  study  is  in  Appendix  D. 

For  direct  comparison  of  the  different  propellants,  the  major  analysis  kept  a  constant  loading 
density  of  approximately  0.85  g/cm^  for  the  gun  simulation.  The  results  are  summarized  in 
Table  8 — each  propellant  is  listed  with  the  three  grain  types,  their  corresponding  web  sizes  for 
448  MPa  maximum  pressure,  and  the  calculated  projectile  exit  velocity  in  each  situation. 
Propellant  types  with  a  superscript  “a”  after  the  name  included  at  least  one  granulation  satisfying 
the  minimum  user  requirements.  M30A1,  JA2,  and  EX-99  each  met  the  velocity  goal  with  one 
or  more  multiperf  granulations,  while  BAMO  was  the  only  propellant  type  to  qualify  with  all  the 
grain  types.  The  rest  of  the  formulations,  alfliough  they  all  had  potential  for  848  m/s  as  shown  in 
the  constant-pressure  calculations,  could  not  reach  that  level  when  simulated  as  a  standard  charge 
consisting  of  single-,  7-,  or  19-perf  grains.  BAMO  exhibited  the  highest  computed  velocity  in 
both  constant-pressure  and  granular  configurations;  JA2  was  the  next-best  computed  velocity  in 
both  simulations. 


Unlike  the  majority  of  propellants  listed  in  Table  8,  EX-99  shows  higher  computed  velocities 
for  the  7-perf  grain  type  than  for  the  19-perf.  This  is  because  of  “two-tiered”  burning  rate  data 
and  their  effect  in  the  gun  model.  Many  IB  programs  (including  IBHVG2)  represent  the 
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Table  8.  Summary  of  Gun  System  Study  Performance  Calculations 


Propellant 

Grain  Geometry 

19-Perf 

7-Perf 

1-Perf 

MIO 

844  m/s 
(2.56  mm) 

841  m/s 
(2.76  mm) 

807  m/s 
(3.82  mm) 

M26E1 

836  m/s 
(3.48  mm) 

835  m/s 
(3.76  mm) 

806  m/s 
(5.20  mm) 

M30AU 

861  m/s 
(3.06  mm) 

856  m/s 
(3.34  mm) 

812  m/s 
(4.72  mm) 

JA2" 

875  m/s 
(4.42  mm) 

867  m/s 
(4.91  mm) 

812  m/s 
(7.20  mm) 

BAMO* 

881  m/s 
(3.13  mm) 

877  m/s 
(3.47  mm) 

859  m/s 
(5.06  mm) 

HELP-42 

807  m/s 
(0.883  mm) 

813  m/s 
(0.980  mm) 

836  m/s 
(1.375  mm) 

EX-99" 

845  m/s 
(2.29  mm) 

858  m/s 
(2.52  mm) 

845  m/s 
(3.57  mm) 

*  Included  at  least  one  granulation  satisfying  die  minunum  user  requirements. 


propellant  burning  rate  equation  in  the  form  i?=aP”,  where  P  is  the  instantaneous  pressure,  n  is 
its  exponent,  and  a  is  the  coefficient  used  to  calculate  R  as  the  rate  of  linear  regression  of  the 
propellant  surface.  Usually  the  values  of  a  and  n  do  not  change  for  the  entire  range  of  pressures 
encountered  during  the  IB  calculation,  as  is  the  case  for  the  first  five  propellants  in  Table  8 
(although  each  propellant  has  its  own  unique  values  for  a  and  n).  The  EX-99  data  were  supplied 
with  experimental  burning  rate  values — ^three  different  pressures  {Pi,  P2,  P3)  and  a  corresponding 
rate  (r;,  r2,  rs)  for  each.  IBHVG2  used  these  values  to  create  an  aP"  representation  between  the 
first  two  pressures,  and  again  between  the  last  two  pressures.  The  solution  values  for  a  and  n 
were  not  the  same  over  the  two  pressure  regions  for  EX-99.  In  effect,  the  IB  code  used  values  ai 
and  ni  for  pressures  between  Pj  and  P2  (and  for  all  pressure  levels  below  Pj),  and  values  aj  and  nj 
for  pressures  between  P2  and  P3  (and  for  all  pressure  levels  above  P3).  The  “two-tier”  rates  can 
produce  effects  such  as  the  higher  exit  velocity  for  the  7-perf  EX-99  charge  than  for  its 
counterpart  19-  and  single-perf  granular  charges  when  the  breakpoint  (P2)  is  in  the  range  of 
pressures  experienced  during  the  IB  cycle. 
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The  exponent  n  for  HELP-42  was  much  higher  than  that  of  the  other  propellants;  its  effect  on 
the  model  was  to  produce  the  highest  projectile  exit  velocity  for  the  single-perf  Help-42 
propellant  and  the  lowest  velocity  for  its  corresponding  19-perf  charge.  This  result  is  an 
interaction  with  the  projectile/bore  resistance  profile  and  the  timing  of  gas  production  by  the 
burning  granular  surface. 

3.7  Propellant  Zoning.  Although  there  was  no  requirement  for  a  lower-velocity  charge  in 
this  study,  a  calculation  of  lighter  charge  weights  for  the  JA2  19-perf  grain  solution  gave  the 
information  in  Table  9.  The  “Charge  Burned”  column  shows  the  fraction  of  main  charge 
consumed  at  the  time  of  projectile  exit — none  of  the  charge  weights  smaller  than  1 1.0  kg  burned 
out.  If  a  low-zone  artillery  application  had  been  a  requirement  with  velocity  at  or  below 
760  m/s,  then  this  propellant  granulation  (web  of  4.43  mm  and  no  changes  in  either  length  or 
perforation  diameter)  would  not  meet  specifications  unless  incomplete  charge  burnout  was 
acceptable  to  the  user. 

3.8  Propellant  Suitability.  The  final  step  in  the  first  round  of  propellant  performance 
testing  could  be  as  simple  as  listing  data  such  as  that  in  Table  8  and  providing  it  to  the  original 
requesting  authority.  Li  this  study,  BAMO  appears  to  be  the  most  suitable  propellant,  since  it 
gives  the  highest  muzzle  energy  and  has  the  flexibility  of  using  any  of  the  three  granulations. 
M30A1  and  JA2  both  provided  solutions  with  the  multiperforated  grain  types  but  fell  short  with 
the  single-perf  geometry.  More  likely,  there  will  be  additional  constraints  put  on  final  selection, 
such  as  the  following: 

•  Propellants  will  be  ranked  according  to  highest  calculated  velocities  within  a  particular 
grain  geometry; 

•  Environmental  and  economic  costs  due  to  manufacturing  or  demilitarization  will 
determine  rankings; 

•  Additional  performance  details,  such  as  computed  chamber  combustion  temperatures, 
will  be  factored  in  with  velocities  to  produce  a  ranking  scheme; 
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Table  9.  Charge  Zoning  Applications  With  19-Perf  JA2  Propellant 


Further  testing,  such  as  supplemental  calculations  with  other  projectile  weights  or 
secondary  propelling  charges,  will  be  indicated; 

Temperature  variations  should  be  considered  if  the  propellant  is  to  be  used  in  climatic 
extremes  of  desert,  tropical,  or  polar  regions; 
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•  If  no  candidate  propellants  survive  the  process,  the  requester  may  ask  for  variations  in 
ingredient  ratios  in  order  to  find  a  viable  charge; 

•  Follow-on  calculations  with  one-,  two-,  or  higher-dimensional  IB  programs  may  be 
indicated  so  that  pressure  waves  or  other  secondary  interior  ballistic  effects  might  be 
discovered  [12, 13]. 

The  comment  concerning  the  choice  of  IB  code  should  be  addressed  before  this  report  could 
be  considered  complete.  IBHVG2  was  used  because  of  its  automatic  input-varying  feature. 
Most  of  the  currently  available  one-  or  multi-dimensional  IB  programs  lack  the  ability  to 
automatically  vary  any  of  the  input  parameters.  In  order  to  complete  several  of  the  process  steps 
(finding  web  size  versus  grain  geometry,  confuting  various  charge  weights  and  noting  amount 
burned  at  time  of  projectile  exit,  etc.),  the  number  of  individual  IB  cycles  to  be  calculated  can  be 
very  large.  The  labor  required  to  modify  input  data,  compute  each  situation,  tabulate  output 
values,  and  analyze  the  results  could  be  extremely  large  if  each  calculation  is  done  separately. 
Therefore,  the  use  of  a  code  with  an  automated  parametric  function  is  crucial,  especially  during 
the  initial  phase  of  propellant  screening.  The  same  type  of  reasoning  can  lead  to  tremendous 
savings  in  time  and  labor  in  virtually  all  steps  of  the  propellant  developmental  process  where 
automated  models  can  be  fitted  with  parameter  variation  [25]. 

3.9  Next  Step.  Other  “tools”  in  the  propellant  production  and  testing  flow  diagram  will 
need  the  information  from  a  performance  analysis,  whether  the  results  are  used  for  re¬ 
formulation  of  a  new  propellant  or  for  additional  filtering  of  existing  choices.  Before  the  user 
elects  an  energetic  material,  he  will  require  more  data  about  manufacturing  costs,  vulnerability, 
shelf  life,  environmental  conpatibility,  and  a  host  of  other  factors  dealing  with  propellants.  The 
initial  performance  comparisons  are  a  small  but  important  part  of  the  overall  formulation  and 
selection  process. 
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4.  Lessons  Learned  and  Tools  Needed 


The  example  test  case  is  a  simplistic  exercise;  no  variations  of  gun  hardware  were 
considered;  no  propellant  reformulation  was  envisioned;  no  multidimensional  IB  analysis  was 
performed.  Yet  it  is  clear  that  some  functions  are  necessary  and  others  are  nonexistent.  To  be 
complete,  the  model-based  performance  process  still  requires  software  tools  to  estimate  density 
and  burning  rates  of  virtual  propellants — semi-empirical  and  first-principle  burning  rate  models 
are  currently  under  development  at  the  Army  Research  Laboratory  [26],  and  models  exist  to 
predict  crystal  densities  within  a  subset  of  crystalline  space  groups  [27].  Automated  search 
methods  for  IB  solutions  are  absolutely  necessary  for  minimization  of  the  labor  involved  in 
computation  and  analysis.  As  the  overall  propellant  creation  and  testing  methods  become 
formalized,  there  will  certainly  be  other  areas  where  software  tools  need  to  be  developed  or 
utilized  for  individual  and  composite  steps  in  the  sequence — estimation  of  processing  costs, 
toxicity  hazards,  mechanical  properties,  disposal  and  demilitarization  requirements,  and  many 
other  activities. 

Along  with  the  computer  programs  is  a  critical  need  for  databases  to  be  assembled  with  the 
parameters  required  for  software  tools.  With  enough  data,  the  missing  tools  for  estimation  of 
propellant  density  and  burning  rates  might  be  created;  vulnerability,  compatibility,  temperature 
sensitivity,  and  a  dozen  other  empirical  methods  could  follow.  At  the  least,  such  gathered  data 
would  make  the  tasks  of  modelers  much  easier  and  faster. 

Clear  and  complete  performance  and  environmental  requirements  might  shorten  the 
performance  modeling  too.  A  standardized  questionnaire  could  alert  the  modeler  to  additional 
factors  in  zoning,  temperature  maximums,  projectile  acceleration  limits,  and  other  stipulations 
not  apparent  during  construction  of  an  initial  set  of  performance  requirements.  Ideally,  a 
database  of  existing  environmental  regulations  and  constraints  should  be  incoiporated  into  the 
process  to  flag  the  user  about  possible  problems  with  manufacture,  movement,  and  disposal  of 
energetic  materials. 
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5.  Summary 


We  have  demonstrated  the  preliminary  steps  in  a  new  formulation  process  for  gun 
propellants  that  can  minimize  life-cycle  costs  by  science-based  design.  The  earhest  steps  in  this 
process  use  modeling  and  simulation  of  the  properties,  processing,  and  performance  of  the 
candidate  in  order  to  screen  potential  candidates,  eliminate  probable  poor  performers,  and 
identify  candidates  that  show  promise.  In  this  exercise,  user  requirements  of  maximum  pressure 
limitation  and  minimum  muzzle  energy  for  a  Navy  5-inch  54-caliber  gun  were  presented.  The 
most  promising  candidate(s)  to  meet  the  user  requirement  from  a  set  of  10  existing  formulations 
representing  single-,  double-  and  triple-based  propellants,  as  well  as  other  novel  energetic 
mixtures,  was  determined.  Since  not  all  performance  and  properties  models  exist,  it  was 
necessary  to  use  limited  empirical  information  to  perform  parts  of  the  exercise.  Down-selection 
of  the  propellants  proceeded  first  through  consideration  of  ingredient  toxicity,  then  calculation  of 
thermodynamic  properties,  estimation  of  theoretical  maximum  energy,  determination  of  actual 
granular  solutions,  and  listing  of  additional  final  considerations.  At  the  end  of  the  exercise,  4  of 
the  10  original  formulations  survived  the  simulation,  and  information  is  available  at  this  point  for 
a  user  to  rank  the  suitability  of  the  propellant  according  to  his  requirements. 

The  number  of  steps  in  a  complete  screening  analysis  would  seemingly  be  endless,  just  as  the 
infinite  amount  of  propellant  variations  to  be  considered  might  appear  when  first  contemplated. 
But  the  hst  of  available  propellant  formulations  is  small  when  compared  to  current  con^uter 
capabilities,  and  a  performance  testing  sequence  with  defined  steps  and  parametric  searching 
methods  can  quickly  determine  which  propellants  can  meet  user  requirements  without  extensive 
firing  range  exercises.  When  compared  to  cost  of  experimental  tests  of  the  same  materials,  the 
amount  of  resources  saved  by  model-based  filtering  is  enormous. 


21 


Intentionally  left  blank. 


22 


6.  References 


1.  Cohen,  W,  S.  “Report  of  the  Quadrennial  Defense  Review.”  U.S.  Department  of  Defense, 
http://www.defenselink.mil/quad.html,  May  1997. 

2.  Miller,  M.  S.,  B.  M.  Rice,  and  R.  J.  Cramer.  “A  New  Approach  to  Propellant  Formulation: 
Minimizing  Life-Cycle  Costs  Through  Science-Based  Design.”  Memorandum  report, 
U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  to  be  published. 

3.  Baer,  P.,  and  J.  Frankie.  “The  Simulation  of  Interior  Ballistic  Performance  of  Guns  by 
Digital  Computer  Program.”  BRLR  1183,  U.S.  Army  Ballistic  Research  Laboratory, 
Aberdeen  Proving  Ground,  December  1962. 

4.  Voltis,  P.  M.  “Digital  Computer  Simulation  of  the  Interior  Ballistic  Process  in  Guns.” 
TR-6615,  Watervliet  Arsenal,  Waverliet,  NY,  1966. 

5.  Baer,  P.  G.,  I.  W.  May,  and  J.  M.  Frankie.  “A  Comparison  of  Several  Predictive 
Approaches  in  Charge  Establishment  for  Large  Caliber  Artillery  Systems.”  11th  JANNAF 
Combustion  Meeting,  Chemical  Propulsion  Information  Agency,  publ.  261,  vol.  1, 
pp.  55-66, 1974. 

6.  Fried,  L.  E.,  W.  M.  Howard,  and  P.  C.  Souers.  “Cheetah  2.0  User’s  Manual.”  Energetic 
Materials  Center,  Lawrence  Livermore  National  Laboratory,  Livermore,  CA,  August  1998. 

7.  Freedman,  E.  “BLAKE— A  Thermodynamics  Code  Based  on  TIGER:  Users’  Guide  to  the 
Revised  Program.”  ARL-CR-422,  U.S.  Army  Research  Laboratory,  Aberdeen  Proving 
Ground,  MD,  July  1998. 

8.  Svehla,  R.  A.,  and  B.  J.  McBride.  “FORTRAN  IV  Computer  Program  for  Calculation  of 
Thermodynamic  and  Transport  Properties  of  Complex  Chemical  Systems.”  NASA 
TD-7056,  NASA  Lewis  Research  Center,  Cleveland,  OH,  1973. 

9.  McBride,  B.  J.,  and  S.  Gordon.  “Computer  Program  for  Calculation  of  Complex  Chemical 
Equilibrium  Compositions  and  Applications,  H:  User’s  Manual  and  Program  Description.” 
NASA  RP-131 1,  NASA  Lewis  Research  Center,  Cleveland,  OH,  1996. 

10.  Anderson,  R.  D.,  and  K.  D.  Fickie.  “IBHVG2— A  User’s  Guide.”  BRL-TR-2829, 
U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  July  1987. 

11.  Oberle,  W.  F.  “Constant  Pressure  Interior  Ballistics  Code  CONPRESS:  Theory  and  User’s 
Manual.”  ARL-TR-199,  U.S.  Army  Research  Laboratory,  Aberdeen  Proving  Ground,  MD, 
September  1993. 


23 


12.  Gough,  P.  S.  “The  XNOVAKTC  Code.”  BRL-CR-627,  U.S.  Amy  Ballistic  Research 
Laboratory,  Aberdeen  Proving  Ground,  MD,  February  1990. 

13.  Nusca,  M.  J.  (Amy  Research  Laboratory),  and  P.  S.  Gough  (Paul  Gough  Associates). 
“Numerical  Model  of  Multiphase  Flows  Applied  to  Solid  Propellant  Combustion  in  Gun 
Systems.”  AIAA  98-3695,  34th  AIAA/ASME/SAE/ASEE  Joint  Propulsion  Conference  & 
Exhibit,  Cleveland,  OH,  13-15  July  1998. 

14.  Mallincktrodt  Baker,  Inc.  “Lead  Carbonate — ^MSDS  Number  L2650.”  Material  Safety  Data 
Sheet,  http://www.jtbaker.com/msds/12650.htm,  December  1996. 

15.  International  Chemical  Safety  Cards.  “2,4-Dinitrotoluene — ^ICSC:  0727.” 

http://hazard.com/msds/mf/cards/file/0727.html,  1 993. 

16.  Freedman,  E.  “Themodynamic  Properties  of  Gun  Propellants.”  Chapter  5  of  Gun 
Propulsion  Technology,  edited  by  L.  Stiefel  and  M.  Summerfield,  vol.  109  of  Progress  in 
Astronautics  and  Aeronautics,  American  Institute  of  Aeronautics  and  Astronautics,  Inc., 
pp.  103-132,  1988. 

17.  Miller,  M.  S.  “Themophysical  Properties  of  Six  Solid  Gun  Propellants.”  ARL-TR-1322, 
U.S.  Amy  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  March  1997. 

18.  Almeyda,  N.  Personal  communication.  Naval  Surface  Warfare  Center,  Indian  Head,  MD, 
November  1998. 

19.  Cramer,  R.  J.  “Advanced  Gun  Propellants.”  CPIA  Ihiblication  675,  vol.  n,  pp.  145-148, 
July  1998. 

20.  Virginia,  F.  J.,  Jr.  “Interior  Ballistics  Firing  Library  of  Closed  Breech,  Single  Combustion 
Chamber  Tank  Guns,  Artillery  Guns,  and  Howitzers  Cannon.”  Report  87-027,  U.S.  Amy 
Research,  Development,  and  Engineering  Center,  March  1987. 

21.  Van  der  Waals,  J.  D.,  Sr.  “Over  de  Continuiteit  van  Den  Gas-en  Vloestoftoestand.”  Ph.D. 
dissertation,  Leiden,  Netherlands,  1873. 

22.  Robbins,  F.  W.,  and  F.  R.  Lyim.  “Analytic  Solutions  to  the  Closed  Bomb.”  BRL-TR-2892, 
U.S.  Amy  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  March  1988. 

23.  Oberle,  W.  F.,  IB,  and  D.  E.  Kooker.  “BRLCB:  A  Closed  Chamber  Data  Analysis  Program 
With  Provisions  for  Deterred  and  Layered  Propellants.”  BRL-TR-3227,  U.S.  Amy  Ballistic 
Research  Laboratory,  Aberdeen  Proving  Ground,  MD,  April  1991. 

24.  Price,  C.,  and  A.  Juhasz.  “A  Versatile  User-Oriented  Closed  Bomb  Data  Reduction 
Program  (CBRED).”  BRL-TR-2018,  U.S.  Amy  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  MD,  September  1977. 


24 


25.  Murphy,  Michael  J.  “GLO — Global  Local  Optimizer  User’s  Manual.”  UCRL-MA-133858, 
Lawrence  Livermore  National  Laboratory,  Livermore,  CA,  April  1999. 

26.  Miller,  M.  S.,  and  W.  R.  Anderson.  “Energetic-Material  Combustion  Modeling  With 
Elementary  Gas-Phase  Reactions:  A  Practical  Approach.”  Journal  of  Propulsion  and 
Power,  to  be  published. 

27.  Holder,  J.  R.,  Z.  Du,  and  H.  L.  Ammon.  “Prediction  of  Possible  Crystal  Structures  for  C-, 
H-,  N-,  0-,  and  F-Containing  Organic  Compounds.”  Journal  of  Computational  Chemistry, 
vol.  14,  no.  4,  p.  422, 1993. 


25 


Intentionally  left  blank. 


26 


Appendix  A: 

BLAKE  Calculations  for  Generic  Ml  Propellant 
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*★**★★*★***★★  BLAKE  Thermodynamic  Equilibrium  &  Gun  Code  ****************** 

*  ARLware  <  Director,  WMRD/ARL,  APG,  MD  21005-5066  * 

*  Attn  :  A.  Kotlar,  AMSRL-WM-BD  >  * 


***  BLAKE  Version  221.50  *** 

Tyger!  Tyger!  Burning  bright/  In  the  forests  of  the  night. 

What  immortal  hand  or  eye/  Dare  frame  thy  fearful  symmetry? 

---William  Blake  (1757-1827) 

16:07:20  6-APR-1999 

>  TITLE,  GENERIC  Ml 

>  COM,NC1315,83.11,DNT,9.77,DBP,4.89,DPA,0.98,H2O,0 .50,ALC,0.75 


«<  The  binary  library  being  used  is  based  on  »> 
«<  SBLAKLYB.LIB  dated  11  August  1998  »> 
«<  »> 
«<  This  binary  library  was  created  on  24 -SEP- 1998  »> 


16:07:20  6-APR-1999  Page  1. 

GENERIC  Ml 

THE  COMPOSITION  IS 


Name 

Pet  Wt 

Pet  Mole 

Delta  H  (J/mol) 

FORMULA 

NC1315 

83.110 

.244 

-6.8890E+08 

C 

H 

0 

N 

6000 

7364 

10271 

2636 

DNT 

9.770 

44.210 

-6.5900E+04 

C 

H 

0 

N 

7 

6 

4 

2 

29 


DBP  4.890  14.480 


DPA  .980  4.773 


H20  .500  22.875 


ALC  .750  13.418 


-8.4260E+05 


1.1190E+05 


-2.8583E+05 


-2.7710E+05 


C  H  0 

16  22  4 

C  N  H 

12  1  11 

H  0 

2  1 

C  H  0 

2  6  1 


The  Elements  and  their  Atom  Percentages 
C  25.597 

H  31.357 

0  34.025 

N  9.022 

Formula  Weight  =  1009.775 

The  Heat  of  Formation  is  -2340.7  J/g  =  -2.364E+06  J/mol 

=  -559.46  cal/g  =  -5.649E+05  cal/mol 


>  GUN, 0.05, 0.05, 0.4 


CONSTITUENT 

CONCENTRATIONS  -  MOLES  PER 

KGM  OF  COMPOUND 

NAME 

1) 

2) 

3) 

4) 

CO 

GAS 

2 .29775E+01 

2.29841E+01 

2.29721E+01 

2.29376E+01 

H2 

GAS 

9.50757E+00 

9.44288E+00 

9 . 35897E+00 

9.24995E+00 

H20 

GAS 

5.98865E+00 

6.02345E+00 

6.05732E+00 

6.09215E+00 

N2 

GAS 

4 .46041E+00 

4 .45214E+00 

4 .44157E+00 

4.42848E+00 

C02 

GAS 

2.36138E+00 

2 .33947E+00 

2.32619E+00 

2.32276E+00 

H 

GAS 

1 . 69823E-02 

1.12191E-02 

8 . 56769E-03 

6.96643E-03 

OH 

GAS 

2 .06596E-03 

1.36908E-03 

1.05132E-03 

8. 63253E-04 

NH3 

GAS 

8 .10015E-03 

1.81443E-02 

3.03683E-02 

4.49315E-02 

HCN 

GAS 

4 .77970E-03 

1.10574E-02 

1.91618E-02 

2.94333E-02 

CH4 

GAS 

1.62997E-03 

7.66132E-03 

2.01363E-02 

4.13460E-02 

0 

GAS 

2.48982E-06 

1.13317E-06 

6.92114E-07 

4.82439E-07 

02 

GAS 

6.37302E-07 

2 . 90323E-07 

1.77820E-07 

1.24651E-07 

NO 

GAS 

9 .00247E-05 

6.15728E-05 

4 . 88600E-05 

4.15520E-05 

CHO 

GAS 

6 . 85851E-04 

1.03823E-03 

1.36422E-03 

1.69659E-03 

CH20 

GAS 

1.45379E-03 

3.21654E-03 

5.33672E-03 

7.86246E-03 

30 


HNCO 

GAS 

4.06787E-04 

9.26268E-04 

1.58314E-03 

2.40638E-03 

NH2 

GAS 

4.64501E-05 

6.85664E-05 

8.76048E-05 

1.05662E-04 

CHS 

GAS 

7.31242E-05 

2.39203E-04 

5.06736E-04 

8.95092E-04 

NH 

GAS 

6.03941E-07 

6.11197E-07 

6.19887E-07 

6.33044E-07 

C2H2 

GAS 

1.25930E-05 

6.41052E-05 

1.83132E-04 

4.10910E-04 

HNO 

GAS 

3.97127E-07 

4.00910E-07 

4.06277E-07 

4.15525E-07 

H02 

GAS 

1.21678E-08 

8.17384E-09 

6.38757E-09 

5.38038E-09 

C2H4 

GAS 

1.04761E-06 

1.10999E-05 

4.93235E-05 

1.52020E-04 

N 

GAS 

1.00572E-07 

6.93893E-08 

5.55090E-08 

4.75971E-08 

CN 

GAS 

5.84721E-07 

9.18693E‘-07 

1.25311E-06 

1.61986E-06 

N20 

GAS 

3.61315E-08 

3.74717E-08 

3.90376E-08 

4.10817E-08 

NCO 

GAS 

2,73096E“07 

4.28349E“07 

5.84300E-07 

7.57305E’07 

HN02 

GAS 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

CH2 

GAS 

2.55410E-07 

5,74641E-07 

9.69875E-07 

1.45750E'06 

C 

GAS 

0 .OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

N02 

GAS 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

FORMAC 

GAS 

5.42035E--04 

1.19880E-03 

1.99208E-03 

2.94845E-03 

C(S) 

SOLID 

O.OOOOOE+00 

O.OOOOOE+00 

0 .OOOOOE+00 

0 . OOOOOE+00 

TOTAL 

GAS  (MOLES/KG)  45.3324 

45.2984 

45.2466 

45.1710 

CONSTITUENT 

CONCENTRATIONS  -  MOLES  PER 

KGM  OF  COMPOUND 

NAME 

5) 

6) 

7) 

8) 

CO 

GAS 

2.28777E+01 

2.27878E+01 

2.26665E+01 

2.25136E+01 

H2 

GAS 

9.11047E+00 

8.93571E+00 

8.72471E+00 

8.47960E+00 

H20 

GAS 

6.12997E+00 

6.17114E+00 

6.21603E+00 

6.26414E+00 

N2 

GAS 

4.41323E+00 

4.39517E+00 

4.37445E+00 

4.35116E+00 

C02 

GAS 

2.33074E+00 

2.35106E+00 

2.38457E+00 

2.43146E+00 

H 

GAS 

5.88089E-03 

5.09719E-03 

4.50813E-03 

4.05113E-03 

OH 

GAS 

7.40195E-04 

6.56138E-04 

5.97905E-04 

5.57619E-04 

NH3 

GAS 

6.18910E-02 

8.11497E-02 

1.02486E-01 

1.25583E-01 

HCN 

GAS 

4.21954E-02 

5.77090E-02 

7.61666E-02 

9.76892E-02 

CH4 

GAS 

7.34602E-02 

1.17984E-01 

1.75340E-01 

2.44728E-01 

0 

GAS 

3.65674E-07 

2.95264E-07 

2.51069E-07 

2.22888E-'07 

02 

GAS 

9.53861E-08 

7.81214E*'08 

6.77253E-08 

6.16155E-08 

NO 

GAS 

3.70115E-05 

3.42002E-05 

3.26063E-05 

3.19280E-05 

CHO 

GAS 

2.05200E-03 

2.44228E-03 

2.87740E-03 

3.36568E-03 

CH20 

GAS 

1.08379E-02 

1,42960E*’02 

1.82605E-02 

2.27449E-02 

HNCO 

GAS 

3.42923E-03 

4.68750E-03 

6.22034E-03 

8.06945E-03 

NH2 

GAS 

1.23644E-04 

1.42049E-04 

1.61195E-04 

1.81242E-04 
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CH3 

GAS 

1.42466E-03 

2.11222E-03 

2.96843E-03 

3 .99649E-03 

NH 

GAS 

6.52591E-07 

6.79944E-07 

7.16152E-07 

7.61697E-07 

C2H2 

GAS 

8 .02648E-04 

1.42623E-03 

2.35908E‘03 

3 .68390E-03 

HNO 

GAS 

4.30369E-07 

4.52266E-07 

4.82529E-07 

5.22183E-07 

H02 

GAS 

4 ,78366E-09 

4.45251E-09 

4.31702E-09 

4 .33842E-09 

C2H4 

GAS 

3 .79141E’’04 

8.17489E-04 

1.57827E-03 

2.78883E-03 

N 

GAS 

4.27529E-08 

3.98361E-08 

3.82821E-08 

3.77544E-08 

CN 

GAS 

2.03988E-06 

2.53223E-06 

3.11645E-06 

3.81221E'06 

N20 

GAS 

4 ,38331E“08 

4.75197E-08 

5.23836E-08 

5.86637E-08 

NCO 

GAS 

9.59587E-07 

1.20374E-06 

1.50417E-06 

l,87724E-06 

HN02 

GAS 

O.OOOOOE+00 

O.OOOOOE+00 

1.80811E-09 

2.04088E-09 

CH2 

GAS 

2,05655E-06 

2.78752E-06 

3.67181E-06 

4 .73006E-06 

C 

GAS 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+OO 

N02 

GAS 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+00 

O.OOOOOE+OO 

FORMAC 

GAS 

4 .10009E-03 

5.48343E-03 

7.14035E-03 

9.11716E-03 

C(S) 

SOLID 

0 . OOOOOE+00 

O.OOOOOE-f-00 

O.OOOOOE+OO 

O.OOOOOE+OO 

TOTAL 

GAS  (MOLES/KG)  45.0694 

44 . 9350 

44.7670 

44.5666 

*  *  SUMMARY  OF  PROPELLANT  THERMO  PROPERTIES  *  * 


Truncated  virial  equation  of  state  with  L-J  6,12  potential  is  being  used 


Rho/L 

Temp 

Press 

Imptus 

Mol  Wt 

Co-Vol 

Frozen 

Balrgy 

U 

PHI 

g/cc 

K 

MPa 

J/g 

Gas 

cc/g 

Gamma 

J/g 

J/g 

1)  .05000 

2436. 

48.79 

918.2 

22.059 

1.181 

1,2630 

3490.9 

-2340.7 

1.0627 

2) .10000 

2439. 

103.9 

918.7 

22.076 

1.158 

1.2639 

3480.9 

-2340.7 

1.1310 

3) .15000 

2443  . 

166.1 

919.0 

22.101 

1.133 

1.2654 

3462.4 

-2340,7 

1.2047 

4)  .20000 

2447. 

236.0 

919.2 

22.138 

1.105 

1.2675 

3436.0 

-2340.7 

1.2838 

5)  .25000 

2454  . 

314.5 

919.5 

22.188 

1,077 

1.2702 

3402.9 

-2340.7 

1.3683 

6)  .30000 

2462. 

402.4 

919.7 

22.254 

1.048 

1.2734 

3364.1 

-2340.7 

1.4583 

7) .35000 

2471. 

500.2 

919,9 

22.338 

1.018 

1.2770 

3320.9 

-2340.7 

1.5536 

8) .40000 

2483. 

608.9 

920.1 

22.438 

.989 

1.2810 

3274 . 3 

-2340.7 

1.6545 

>  QUIT 

Run  time  *  1.15  seconds 
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Appendix  B: 

IBHVG2  Constant-Pressure  Calculations 
Using  M26E1  Propellant 
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Intentionally  left  blank. 
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ERRTOL-  1.1920929E-07 

1 


IBHVG2.505 


DATE 


TIME 


0  CARD 

1  -'> 

$COMM 

CARD 

2  --> 

GREEN  GUN  TEST  CASE  -  INDIAN  HEAD  SYSTEM 

CARD 

3  --> 

$GUN 

CARD 

4  -“> 

TRAV  =  6.840  CHAM  -  0.01506  $  919  CUBIC  INCHES 

CARD 

5  --> 

GRVE  =  0.12852  LAND  -  0.127  TWST  -  24  G/L  =  1.66  $ 

ESTIMATES 

CARD 

6  --> 

$  CRTS  -  6 

CARD 

7  --> 

$  DIST  ■=  0.0,  0.04455,  0.14732,  0.8260,  0.8913,  1.0592 

CARD 

8 

$  DIAM  =  0.12965,  0.13655,  0.13929,  0.132588,  0.12852, 

0.12852 

CARD 

9  --> 

$PROJ 

CARD 

10  --> 

PRWT  =  49.895 

CARD 

11  --> 

$RESI 

CARD 

12  *-> 

NPTS  =  5 

CARD 

13  --> 

TRAV  =  0.0,  0.00508,  0.04318,  0.1397,  7.62 

CARD 

14  --> 

PRES  =  3.4474,  0.6894757,  0.6894757,  8.2737,  8.2737 

CARD 

15  --> 

$INFO 

CARD 

16  --> 

POPT  =  1,1, 1,0 

CARD 

17  --> 

RUN  =  ’GREEN  GUN  TEST  CASE  -  M26E1’ 

CARD 

18  --> 

$  GRAD  -  3 

CARD 

19  --> 

CONP  -  2  PRES  =  448 

CARD 

20  --> 

$COMM  PRIM 

CARD 

21  --> 

NAME  =  ’PRIMER'  CHWT  =  0.12637  FORC  *  146140 

CARD 

22  --> 

GAMA  =  1.2015  COV  =  0.0010838  TEMP  “  3377 

CARD 

23  --> 

$PROP 

CARD 

24 

NAME  =  'M26E1'  CHWT  =  11.567  FORC  =  1085031 

CARD 

25 

GAMA  *=  1.2384  COV  =  0.0010383  TEMP  *=  3132 

CARD 

26  --> 

FORM  -  '7P' 

CARD 

27  --> 

LEN  *  0.018288  PD  =  0.0004572  WEB  =  0.003467  RHO  «= 

1611 

CARD 

28  --> 

ALPH  =  0.7468  BETA  =  0.0033778 

CARD 

29  --> 

$COMM  PARA 

CARD 

30 

VARY  =  ’CHWT’  DECK  *=  'PROP'  NTH  1 

CARD 

31  --> 

FROM  *  11.5  BY  =  0.1  TO  =  15.01 

CARD 

32  --> 

$END 

IGREEN 

GUN  TEST  CASE  -  M26E1  IBHVG2.505 

DATE 

GUN  TUBE 


TYPE: 

(M) :  6.84000 

CHAMBER  VOLUME  (M3): 

0.01506 

GROOVE  DIAMETER  (M)  :  0.12852 

GROOVE/LAND  RATIO  (-):  1.660 

LAND  DIAMETER  (M)  : 

TIME 


TRAVEL 

0.12700 
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TWIST  (CALS/TURN):  24.0  BORE  AREA  (M2):  0.01286 

OPTION:  1 

SHELL  THICKNESS  (M)  :  0.000102  SHELL  CP  (J/KG-K) :  460.3163 

DENSITY  (KG/M3):  7861.0918 

INITIAL  SHELL  TEMP  (K)  :  293.  AIR  HO  (W/M**2-K) :  11,3482 


“  PROJECTILE  - 


TYPE: 

PREDICTOR  OPTION:  0 


TOTAL  WEIGHT  (KG):  49.895 


-  RESISTANCE  - 


AIR  RESISTANCE  OPTION:  1  WALL  HEATING  FRACTION:  0.000 

RESISTIVE  PRESSURE  MULT  INDEX:  3  RESISTIVE  FACTOR  1,000 

TABLE  LENGTH :  5 


I 

TRAVEL  (M) 

PRESSURE  (MPA) 

I  TRAVEL  (M) 

PRESSURE  (MPA) 

TRAVEL 

(M) /  PRESSURE 

(MPA) 

1 

0.000 

3.447 

3  0.043 

0.689 

7.620 

8.274 

2 

0.005 

0.689 

4  0.140 

8.274 

-  GENERAL  - 


MAX  TIME  STEP  (S): 

0.000100 

PRINT  STEP  (S) : 

0.000000 

RELATIVE  ERROR  (-): 

0.00200 

PRINT  OPTIONS: 

1110  11 

STORE  OPTION: 

0 

PRESSURE  OPTION: 

2 

GRADIENT  MODEL:  LAGRANGIAN 

-  RECOIL  - 


RECOIL  OPTION: 

WEIGHT  ( KG ) :  0 . 


0  TYPE : 


HEAT -LOSS 

SHELL 


WEIGHT 


FRICTION 

I 

5 


MAX 

CONSTANT - 


RECOILING 
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CONSTANT -PRESSURE  RUN 


MAINTAIN  PRESSURE  (MPA) ; 
VARYING:  PROPELLANT  SURFACE  AREA 
IGREEN  GUN  TEST  CASE  -  M26E1 


WITHIN  (MPA) : 


IBHVG2.505 


-  CHARGE  1  - 


TYPE:  M26E1 

(KG):  11.5670 

EROSIVE  COEFF  (-): 

AT  (S):  O.OOOOOE+00 


1.0000  7P 


0,000000  CHARGE  IGN  CODE: 


CHARGE  IGN 


PROPERTIES  AT  LAYER  BOUNDARIES  OF  LAT  SURFACES 


AT  DEPTH  (M) : 

ADJACENT  LAYER  WT  %: 
DENSITY  (KG/M3): 

GAMMA  ( - ) ; 

FORCE  (J/KG): 

COVOLUME  (M3/KG): 

FLAME  TEMP  (K) : 

BURNING  RATE  EXPS : 

BURNING  RATE  COEFFS : 
IGREEN  GUN  TEST  CASE  -  M2 6 El 


0.00000 
.  100.000 
1611.000 
1.2384 
1085031. 
-  1.0383E-03 

3132.0 
0.7468 
—  3.3778E-03 

IBHVG2.505 


TIME 

TRAV 

VEL 

ACC 

BREECH 

MEAN 

BASE 

MEAN 

FRAC 

SURFACE 

BURNING 

(MS) 

(M) 

(M/S) 

(G) 

PRESS 

PRESS 

PRESS 

TEMP 

BURN 

AREA 

RATE 

(MPA) 

(MPA) 

(MPA) 

(K) 

1 

(M**2) 

(M/S) 

0,000 

0,000 

0.00 

10469. 

448.000 

432.608 

401.823 

3132. 

0.233 

0.223 

0.314 

0.100 

0.001 

10.27 

10475. 

448.000 

432.598 

401.794 

3131, 

0.233 

0.056 

0.314 

RREL  RESISTANCE  OVERCOME  -  PROJECTILE  MOVING 

0.200 

0.002 

20.55 

10495. 

448.000 

432.569 

401.708 

3128. 

0.234 

0.169 

0.314 

0.300 

0.005 

30.86 

10527. 

448.000 

432.521 

401.563 

3124. 

0.235 

0.283 

0.314 

0.400 

0.008 

41.19 

10533. 

448.000 

432.513 

401.538 

3118. 

0.237 

0.402 

0.314 

0.500 

0.013 

51.53 

10533. 

448.000 

432.513 

401.539 

3110. 

0.239 

0.519 

0.314 

0,600 

0.019 

61.86 

10533. 

448.000 

432.513 

401.539 

3101. 

0.242 

0.634 

0.314 

0.700 

0,025 

72.19 

10533. 

448.000 

432.513 

401.540 

3091. 

0.245 

0.749 

0.314 

0,800 

0.033 

82.52 

10533. 

448.000 

432.514 

401.540 

3079. 

0.249 

0.864 

0.314 

0,900 

0.042 

92.85 

10533. 

448.000 

432.514 

401.541 

3066. 

0.253 

0.979 

0.314 

1,000 

0.052 

103.17 

10517. 

448.000 

432.537 

401.610 

3052. 

0.258 

1.097 

0.314 
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1.212  0.314 
1.326  0.314 
1.440  0.314 


1.100  0.062  113.48  10497.  448.000  432.566  401.699  3038.  0.263 

1.200  0.074  123.77  10475.  448.000  432.599  401.796  3022.  0.269 

1.300  0.087  134.03  10451.  448.000  432.634  401.902  3006.  0.275 

1.400  0.101  144,27  10425.  448.000  432.672  402,017  2990.  0.282  1.554  0.314 

1.500  0,116  154.48  10397.  448.000  432.713  402.139  2973.  0.290  1.667  0.314 

1.600  0.132  164.66  10367.  448.000  432.757  402.270  2957.  0.297  1.779  0.314 

1.700  0,149  174.82  10353.  448.000  432.778  402.334  2940.  0.306  1.888  0.314 

1.800  0.167  184.97  10353.  448.000  432.778  402.335  2923.  0.314  1.997  0.314 

1.900  0.186  195.13  10352.  448.000  432.779  402.336  2906.  0.324  2.107  0.314 

2.000  0.206  205.28  10352.  448.000  432.779  402.337  2890.  0.333  2.218  0.314 

2.100  0.227  215.43  10352.  448.000  432.779  402.337  2874.  0.343  2.328  0.314 

2.200  0.249  225.59  10352.  448.000  432.779  402.338  2858.  0.354  2.439  0.314 

2.300  0,272  235.74  10352.  448.000  432.780  402.339  2842.  0.365  2.549  0.314 

2.400  0.296  245.89  10351.  448.000  432.780  402.340  2827.  0.377  2.659  0.314 

2.500  0.321  256.05  10351.  448.000  432.780  402.341  2812.  0.389  2.768  0.314 

2.600  0.347  266.20  10351.  448.000  432.781  402.342  2798.  0.402  2.878  0.314 

2.700  0.374  276.35  10351.  448.000  432.781  402.343  2784.  0.415  2.987  0.314 

2.800  0.403  286.50  10351.  448.000  432.781  402,344  2771.  0.428  3.096  0,314 

2.900  0.432  296.65  10350.  448.000  432.782  402.345  2758.  0,442  3.206  0.314 

3.000  0.462  306.80  10350.  448.000  432.782  402.346  2745.  0,457  3.315  0.314 

3.100  0.493  316.96  10350.  448.000  432.782  402.347  2733.  0,472  3,424  0.314 

3.200  0.525  327.11  10350.  448.000  432.783  402.348  2722.  0.487  3.533  0.314 

3.300  0.559  337.26  10349.  448.000  432.783  402.349  2711.  0.503  3.642  0.314 

3.400  0,593  347.41  10349.  448.000  432.783  402.350  2700.  0,520  3.750  0,314 

3.500  0.628  357.56  10349.  448.000  432.784  402.352  2690.  0.537  3.859  0.314 

3.600  0.664  367.71  10349.  448.000  432.784  402.353  2680.  0.554  3.968  0.314 

3.700  0.702  377.86  10348.  448.000  432.785  402.354  2671.  0.572  4.077  0.314 

3.800  0.740  388.00  10348.  448.000  432.785  402.355  2662.  0.590  4.186  0.314 

3.900  0.779  398.15  10348.  448.000  432.785  402.356  2653.  0.609  4.295  0.314 

4.000  0.820  408.30  10347.  448.000  432.786  402.358  2645.  0,628  4.404  0,314 

4.100  0.861  418.45  10347.  448.000  432.786  402.359  2637.  0.648  4.513  0.314 

4.200  0.903  428.60  10347.  448.000  432.787  402.360  2629.  0.668  4.622  0.314 

4.300  0.947  438.74  10346.  448.000  432.787  402.362  2622.  0.689  4.731  0.314 

4.400  0.991  448.89  10346.  448.000  432.788  402.363  2615.  0.710  4.841  0.314 

4.500  1.036  459.04  10346.  448.000  432.788  402.365  2608.  0.732  4.950  0.314 

4.600  1.083  469.18  10345.  448.000  432.789  402.366  2602.  0.754  5.059  0.314 

4.700  1.130  479,33  10345.  448.000  432.789  402.368  2595.  0.777  5.169  0.314 

4.800  1.179  489.48  10345.  448.000  432.790  402.369  2589.  0.800  5.278  0.314 

4.900  1.228  499.62  10344.  448.000  432.790  402.371  2584.  0.823  5.388  0.314 

5.000  1.279  509.77  10344.  448.000  432.791  402.372  2578.  0.847  5.497  0.314 

5.100  1.330  519.91  10344.  448.000  432.791  402.374  2573.  0.872  5.607  0.314 

5.200  1.383  530.06  10343.  448.000  432.792  402.376  2568.  0.897  5.717  0.314 

5.300  1.436  540.20  10343.  448.000  432.792  402.377  2563.  0.922  5.827  0.314 

IGREEN  GUN  TEST  CASE  -  M26E1  IBHVG2.505  DATE  T 

TIME  TRAV  VEL  ACC  BREECH  MEAN  BASE  MEAN  FRAC  SURFACE  BURNING 
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(MS) 

(M) 

(M/S) 

(G) 

PRESS 

PRESS 

PRESS 

TEMP 

BURN 

AREA 

RATE 

(MPA) 

(MPA) 

(MPA) 

(K) 

1 

(M**2) 

(M/S) 

1  5.400 

1.491 

550.34 

10343, 

448.000 

432.793 

402.379 

2558. 

0.948 

5.937 

0.314 

5.500 

1.546 

560.49 

10342. 

448.000 

432.794 

402.381 

2554. 

0.975 

6.047 

0.314 

5.593 

1.599 

569.96 

10342. 

448,000 

432.794 

402.382 

2550. 

1,000 

0.000 

0.314 

PROPELLANT  1  BURNED  OUT 

5.600 

1.603 

570.65 

10314. 

446.810 

431.645 

401.316 

2548. 

1.000 

0.000 

0.000 

5.700 

1.660 

580.57 

9915. 

429.871 

415.293 

386.138 

2528, 

1.000 

0.000 

0.000 

5.800 

1.719 

590.10 

9536. 

413.800 

399.779 

371.738 

2509. 

1.000 

0,000 

0.000 

5.900 

1.778 

599.27 

9176. 

398,549 

385.057 

358.073 

2490. 

1.000 

0.000 

0.000 

6.000 

1.839 

608.10 

8835. 

384.073 

371.083 

345,103 

2472. 

1.000 

0.000 

0.000 

6.100 

1.900 

616.61 

8511. 

370.330 

357.816 

332.788 

2454. 

1.000 

0.000 

0.000 

6.200 

1.962 

624,80 

8203. 

357.277 

345.215 

321.093 

2436. 

1.000 

0.000 

0.000 

6.300 

2.025 

632.70 

7911. 

344.875 

333.243 

309.981 

2419. 

1.000 

0.000 

0.000 

6.400 

2.088 

640.32 

7633. 

333.087 

321.864 

299.418 

2402. 

1.000 

0.000 

0.000 

6.500 

2.153 

647.68 

7369. 

321.877 

311.043 

289.374 

2385. 

1.000 

0.000 

0.000 

6.600 

2.218 

654.78 

7117, 

311.212 

300.748 

279.819 

2369. 

1.000 

0.000 

0.000 

6.700 

2.284 

661.64 

6878. 

301.060 

290,948 

270.723 

2354. 

1.000 

0.000 

0.000 

6.800 

2.350 

668.27 

6650. 

291.393 

281.615 

262.060 

2338. 

1.000 

0.000 

0.000 

6.900 

2.417 

674.69 

6433. 

282.181 

272.723 

253.807 

2323. 

1.000 

0.000 

0.000 

7.000 

2.485 

680.89 

6225. 

273.398 

264.245 

245.938 

2308. 

1.000 

0.000 

0.000 

7.100 

2.554 

686.90 

6028. 

265.021 

256.158 

238.432 

2294. 

1.000 

0.000 

0.000 

7.200 

2.623 

692.72 

5839. 

257.026 

248.440 

231.268 

2280. 

1.000 

0.000 

0.000 

7.300 

2.692 

698.35 

5659. 

249.391 

241.070 

224.428 

2266. 

1.000 

0.000 

0.000 

7.400 

2.762 

703.82 

5487. 

242.096 

234.028 

217.892 

2252. 

1.000 

0.000 

0.000 

7.500 

2.833 

709.12 

5323. 

235.123 

227,297 

211.644 

2239. 

1.000 

0.000 

0.000 

7.600 

2.904 

714.26 

5165. 

228,453 

220.858 

205.669 

2226. 

1.000 

0.000 

0,000 

7.700 

2.976 

719.25 

5015. 

222.070 

214.697 

199.950 

2213. 

1.000 

0.000 

0.000 

7.800 

3.048 

724.10 

4871. 

215.959 

208.797 

194.474 

2201. 

1.000 

0.000 

0.000 

7.900 

3.121 

728.81 

4733. 

210.103 

203.145 

189.228 

2189. 

1.000 

0.000 

0.000 

8.000 

3.194 

733.38 

4600. 

204.491 

197.727 

184.199 

2177. 

1.000 

0.000 

0.000 

8.100 

3.267 

737.83 

4473. 

199.108 

192.531 

179.376 

2165. 

1.000 

0.000 

0.000 

8.200 

3.341 

742.16 

4351. 

193.943 

187.545 

174.749 

2153. 

1.000 

0.000 

0.000 

8.300 

3.416 

746.37 

4234. 

188.985 

182.759 

170.306 

2142. 

1.000 

0.000 

0.000 

8.400 

3.491 

750.47 

4122. 

184.222 

178.161 

166.039 

2131. 

1.000 

0.000 

0.000 

8.500 

3.566 

754.45 

4014, 

179.645 

173.743 

161.939 

2120. 

1.000 

0.000 

0.000 

8.600 

3.641 

758.34 

3910. 

175.245 

169.495 

157.996 

2109. 

1.000 

0.000 

0.000 

8.700 

3.717 

762.13 

3810. 

171.012 

165.409 

154.204 

2099. 

1.000 

0.000 

0.000 

8.800 

3.794 

765.81 

3714. 

166.938 

161.477 

150.554 

2089. 

1.000 

0.000 

0.000 

8.900 

3.871 

769.41 

3622. 

163.016 

157.691 

147.040 

2078. 

1.000 

0.000 

0.000 

9.000 

3.948 

772.92 

3533. 

159,237 

154.043 

143.655 

2069. 

1.000 

0.000 

0,000 

9.100 

4.025 

776.34 

3447. 

155.596 

150.529 

140.393 

2059. 

1.000 

0.000 

0.000 

9.200 

4.103 

779.68 

3364. 

152.086 

147.140 

137.248 

2049. 

1.000 

0.000 

0.000 

9.300 

4,181 

782.94 

3284. 

148,700 

143.871 

134.215 

2040. 

1.000 

0.000 

0.000 

9.400 

4.260 

786.12 

3207. 

145.432 

140.717 

131.287 

2030. 

1.000 

0.000 

0.000 

9.500 

4.338 

789.23 

3132. 

142.278 

137,672 

128.461 

2021. 

1.000 

0.000 

0.000 
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9.600 

4.417 

792.27 

3060. 

139.231 

134.732 

125.732 

2012. 

1.000 

0.000 

0.000 

9.700 

4.497 

795.23 

2991. 

136.288 

131.890 

123.095 

2003. 

1.000 

0.000 

0.000 

9.800 

4.576 

798.13 

2924. 

133.443 

129.144 

120.546 

1995. 

1.000 

0.000 

0.000 

9.900 

4.656 

800.97 

2859. 

130.692 

126.489 

118.082 

1986. 

1.000 

0.000 

0.000 

10.000 

4.737 

803.74 

2796. 

128.032 

123.921 

115.698 

1978. 

1.000 

0.000 

0.000 

10.100 

4.817 

806.45 

2735. 

125.457 

121.435 

113.392 

1969. 

1.000 

0.000 

0.000 

10.200 

4.898 

809.11 

2676. 

122.964 

119.029 

111.159 

1961. 

1.000 

0.000 

0.000 

10.300 

4.979 

811.70 

2619. 

120.551 

116.699 

108.996 

1953. 

1.000 

0.000 

0.000 

10.400 

5.060 

814.24 

2564. 

118.212 

114.442 

106.902 

1945. 

1.000 

0.000 

0.000 

10.500 

5.142 

816.73 

2511. 

115.947 

112.255 

104.872 

1937. 

1.000 

0.000 

0.000 

10.600 

5.224 

819.17 

2459. 

113.750 

110.135 

102.904 

1930. 

1.000 

0.000 

0.000 

IGREEN  GUN  TEST  CASE  -  M26E1  IBHVG2.505  DATE  TIME 


TIME 

TRAV 

VEL 

ACC  BREECH 

MEAN 

BASE 

MEAN 

FRAC 

SURFACE 

BURNING 

(MS) 

(M) 

(M/S) 

(G)  PRESS 

PRESS 

PRESS 

TEMP 

BURN 

AREA 

RATE 

(MPA) 

(MPA) 

(MPA) 

(K) 

1 

(M**2) 

(M/S) 

0  10.700 

5.306 

821.56 

2409.  111.620 

108.079 

100.996 

1922. 

1.000 

0.000 

0.000 

10.800 

5.388 

823.89 

2360.  109.554 

106.084 

99.145 

1915. 

1.000 

0,000 

0.000 

10.900 

5.470 

826.18 

2312.  107.549 

104.149 

97.349 

1907. 

1.000 

0.000 

0.000 

11.000 

5.553 

828.43 

2267.  105.603 

102.270 

95.605 

1900. 

1.000 

0.000 

0.000 

11.100 

5.636 

830.63 

2222.  103.713 

100.446 

93.912 

1893. 

1.000 

0.000 

0.000 

11.200 

5.719 

832.79 

2179.  101.878 

98.674 

92.268 

1886. 

1.000 

0.000 

0.000 

11.300 

5.803 

834.90 

2136.  100.094 

96.953 

90.670 

1879. 

1.000 

0.000 

0.000 

11.400 

5.886 

836.98 

2096.  98.361 

95.280 

89.118 

1872. 

1.000 

0.000 

0.000 

11.500 

5.970 

839.01 

2056.  96.676 

93.654 

87.609 

1865. 

1.000 

0.000 

0.000 

11.600 

6.054 

841.01 

2017.  95.038 

92.072 

86.141 

1859. 

1.000 

0 . 000 

0.000 

11.700 

6.138 

842.97 

1979.  93.445 

90.534 

84.713 

1852. 

1.000 

0.000 

0.000 

11.800 

6.223 

844.89 

1943.  91.894 

89.037 

83.324 

1846. 

1.000 

0.000 

0.000 

11.900 

6.307 

846.78 

1907.  90.385 

87.581 

81.973 

1839. 

1.000 

0.000 

0.000 

12.000 

6.392 

848.63 

1872.  88.916 

86.163 

80.657 

1833. 

1.000 

0.000 

0.000 

12.100 

6.477 

850.45 

1839.  87.486 

84.782 

79.375 

1827. 

1-000 

0.000 

0.000 

12.200 

6.562 

852.24 

1806.  86.092 

83.437 

78.127 

1820. 

1.000 

0.000 

0.000 

12.300 

6.647 

854.00 

1774.  84.735 

82.127 

76.911 

1814. 

1.000 

0.000 

0.000 

12.400 

6.733 

855.72 

1742.  83.412 

80.850 

75.726 

1808. 

1.000 

0.000 

0.000 

12.500 

6.819 

857.41 

1712.  82.123 

79.606 

74.572 

1802. 

1.000 

0.000 

0.000 

12.525 

6.840 

857.83 

1704.  81.806 

79.300 

74.287 

1801. 

1.000 

0.000 

0.000 

PROJECTILE  EXIT 

IGREEN  GUN  TEST  CASE  -  M26E1  IBHVG2.505  DATE  TIME 


CONDITIONS  AT: 

PMAX 

MUZZLE 

TIME  (MS): 

5.593 

12.525 

TRAVEL  (M) : 

1.5988 

6.8400 

VELOCITY  (M/S) 

569.96 

857.83 

40 


ACCELERATION  (G) :  10342. 
BREECH  PRESS  (MPA):  447.9999 
MEAN  PRESS  (MPA):  432.7940 
BASE  PRESS  (MPA):  402.3822 
MEAN  TEMP  (K) :  2550. 
Z  CHARGE  1  (-) :  1.000 


ENERGY  BALANCE  SUMMARY 

TOTAL  CHEMICAL: 

(1)  INTERNAL  GAS: 

(2)  WORK  AND  LOSSES: 

(A)  PROJECTILE  KINETIC: 

(B)  GAS  KINETIC: 

(C)  PROJECTILE  ROTATIONAL: 

(D)  FRICTIONAL  WORK  TO  TUBE: 

(E)  OTHER  FRICTIONAL  WORK: 

(F)  WORK  DONE  AGAINST  AIR: 

(G)  HEAT  CONVECTED  TO  BORE: 

(H)  RECOIL  ENERGY: 


1704. 
81.8058 
79.2997 
74.2875 ■ 
1801. 
1.000 


JOULE 

% 

52645464. 

100.00 

30269280. 

57.50 

22376184. 

42.50 

18358272. 

34.87 

1418647. 

2.69 

157282. 

0.30 

0. 

0.00 

718839. 

1.37 

69870. 

0.13 

1653274. 

3.14 

0. 

0.00 

LOADING  DENSITY  (KG/M3) :  768.061 
CHARGE  WT/PROJECTILE  WT:  0.232 
PIEZOMETRIC  EFFICIENCY:  0.466 
EXPANSION  RATIO;  6.840 


41 


Intentionally  left  blank. 
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Appendix  C: 

IBHVG2  Performance  Calculations 
for  M30A1  Propellant 
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Intentionally  left  blank. 
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IGREEN  GUN  TEST  CASE  -  M30A1  19P  PERFORMANCE  IBHVG2.505  DATE  TIME 


TYPE:  CHAMBER  VOLUME  (M3):  0.01506  TRAVEL  (M)  ;  6.84000 
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GENERAL 


MAX  TIME  STEP  (S) :  0.000100  PRINT  STEP  (S):  0.000500  MAX  RELATIVE  ERROR  (“):  0.00200 
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PROPERTIES  AT  LAYER  BOUNDARIES  OP  PERP  SURFACES  PROPERTIES  AT  LAYER  BOUNDARIES  OP  END  SURFACES 
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(A)  PROJECTILE  KINETIC;  18820410.  32.87 


(B)  GAS  KINETIC:  1606417.  2.81 
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4760E“02  254.846  674.44  80.294  0.837  6.840 
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TYPE:  CHAMBER  VOLUME  (M3):  0.01506  TRAVEL  (M) :  6.84000 
GROOVE  DIAMETER  (M) :  0.12852  LAND  DIAMETER  (M) :  0.12700  GROOVE/LAND  RATIO  (-):  1.490 
TWIST  (CALS/TURN) :  25.0  BORE  AREA  (M2);  0.01285  HEAT-LOSS  OPTION:  1 
SHELL  THICKNESS  (M) :  0.000102  SHELL  CP  ( J/KG-K) :  460.3163  SHELL  DENSITY  (KG/M3) :  7861.0918 


INITIAL  SHELL  TEMP  (K):  293.  AIR  HO  (W/M**2-K):  11.3482 
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RECOIL  OPTION:  0  TYPE:  RECOILING  WEIGHT  (KG): 
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BURNING  RATE  EXPS:  .  .  .  0.8063  .  .  .  0.8063 


PROPERTIES  AT  LAYER  BOUNDARIES  OF  LAT  SURFACES 


CN  ro  63  m  00 
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0.3280E-02  780.073  941.68  75.676  1.000  1.333 
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4460E-02  485.075  838.59  87.615  1.000  4.965 
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Appendix  D: 

Performance  Calculations  Varying  Charge  Mass 
of  M30A1  Propellant 
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TWIST  (CALS/TURN):  25.0  BORE  AREA  (M2):  0.01285  HEAT-LOSS  OPTION:  1 
SHELL  THICKNESS  (M)  :  0.000102  SHELL  CP  (J/KG-K) :  460.3163  SHELL  DENSITY  (KG/M3):  7861.0918 
INITIAL  SHELL  TEMP  (K) :  293.  AIR  HO  (W/M**2-K):  11.3482 
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RECOIL  OPTION:  0  TYPE:  RECOILING  WEIGHT  (KG) 


TYPE;  PRIMER  GAMMA  (-):  1.2015  FORCE  (J/KG) ;  146140. 
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BURNING  RATE  COEFFS ;  .  .  .  1.9684E-03  -  -  -  1.9684E-03 


ADJACENT  LAYER  WT%:  .  .  .  100.000 
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